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We explore the likelihood that early remains of Earth, Mars, and 
Venus have been preserved on the Moon in high enough concentra- 
tions to motivate a search mission. During the Late Heavy Bom- 
bardment, the inner planets experienced frequent large impacts. 
Material ejected by these impacts near the escape velocity would 
have had the potential to land and be preserved on the surface of 
the Moon. Such ejecta could yield information on the geochemical 
and biological state of early Earth, Mars, and Venus. To determine 
whether the Moon has preserved enough ejecta to justify a search 
mission, we calculate the amount of terran material incident on 
the Moon over its history by considering the distribution of ejecta 
launched from the Earth by large impacts. In addition, we make 
analogous estimates for Mars and Venus. We find, for a well-mixed 
regolith, that the median surface abundance of terran material is 
roughly 7 ppm, corresponding to a mass of approximately 20,000 kg 
of terran material over a 10 x 10-square-km area. Over the same 
area, the amount of material transferred from Venus is 1-30 kg 
and material from Mars as much as 180 kg. Given that the amount 
of terran material is substantial, we estimate the fraction of this 
material surviving impact with intact geochemical and biological 
tracers. © 2002 Elsevier Science (USA) 
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I. INTRODUCTION 


The frequency of both lunar and martian meteorites on the 
Earth indicates that the transfer of planetary material is common 
in the Solar System. Vigorous hydrologic or tectonic cycles, past 
or present, prevent most nearby planetary bodies from serving 
as long-term repositories of this material. The Moon is an im- 
portant exception, however. Strategically located within the in- 
ner Solar System, the Moon has theoretically collected material 
from all of the terrestrial planets since its formation. Lacking an 
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atmosphere and widespread, long-lasting volcanism, the Moon 
has potentially preserved meteorites from Mercury through the 
asteroid belt. While the lack of an atmosphere prevents a soft 
landing on the lunar surface, its low gravity means particles 
with small velocities with respect to the Moon will experience 
relatively low impact velocities. Moreover, unlike on other ter- 
restrial planets, martian, venusian, and terran meteorites, blasted 
off their respective planets 3.9 Gyr during the Late Heavy Bom- 
bardment, should still exist on the surface of the Moon. Such 
meteorites are likely to contain uniquely preserved remains of 
these planets that are not available elsewhere in the Solar Sys- 
tem. In particular, terran meteorites on the Moon may provide a 
substantive geological record for ancient Earth, corresponding 
to or predating the period for which the earliest evidence for life 
exists. 

In considering the argument in favor of searching the Moon for 
terran meteorites, the SNC meteorites represent obvious analogs. 
A considerable amount of martian geoscience is built upon lit- 
tle more than a dozen samples. From them have been inferred 
important constituents of the atmosphere, mantle, and core, the 
extent of interaction between the martian hydrosphere and litho- 
sphere, constraints on martian water abundance, and the nature 
of a martian carbon sink, while the lightly shocked condition of 
the meteorites has stimulated developments in the understanding 
of impact physics (McSween 1994). Terran meteorites have the 
potential to provide similar information, extending and broaden- 
ing Earth’s geologic record for a time period that has otherwise 
left little or no physical evidence. The rocks’ elemental compo- 
sition and mineralogy (in particular, hydration) could be used 
to constrain characteristics of the early crust and mantle, the 
global oxidation state, the extent of planetary differentiation, 
and the availability of water. Volatile inclusions sampling no- 
ble gases, carbon dioxide, and molecular nitrogen could clarify 
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atmospheric origin and evolution and, along with the meteorite 
mineralogy, could provide substantive constraints on early at- 
mospheric concentrations (Bogard and Johnson 1983). Direct 
measurement of the timing, extent, and planetary effects of the 
heavy bombardment by careful dating of terran meteorites is also 
possible and would perhaps be the most robust and significant 
scientific reward of this project. 

In addition to the scientific benefits listed above, which alone 
justify searching the Moon for terran meteorites, a fraction of 
Earth-derived material on the Moon could contain geochemical 
and biological information, in the form of isotopic signatures, 
organic carbon, molecular fossils, biominerals or even, theoret- 
ically, microbial fossils. Again, analogy to SNC meteorites is 
instructive. Questionable interpretations of structures within the 
martian meteorite ALH84001 as microbial fossils (McKay et al. 
1996) and, more recently, evidence supporting deposition there 
of magnetite by biomineralization (Thomas-Keprta et al. 2001, 
Friedmann et al. 2001) have been used to argue that this mete- 
orite contains vestiges of ancient martian life. While such inter- 
pretations remain highly controversial, they support the general 
principle that terran meteorites should be examined for poten- 
tially novel evidence concerning early Earth life. Such evidence 
could substantiate or extend a contested fossil record that begins 
3.5 Gyr (Awramik 1982, Schopf and Packer 1987, Buick 1990, 
Brasier et al. 2002) and geochemical evidence from even earlier 
periods, between 3.7 and 3.85 Gyr (Schidlowski 1988, Mojzsis 
et al. 1996, Rosing 1999; but see Moorbath 2001 and Fedo and 
Whitehouse 2002). 

In addition, the Moon may preserve material not only from 
Earth, but also from Venus and the asteroid belt. The only at- 
tainable record of Venus’ early surface geology, otherwise catas- 
trophically erased 700 myr ago (McKinnon etal. 1997), is proba- 
bly on the Moon. Similarly, a record of the type, characteristics, 
and origins of the heavy bombardment impactors themselves 
may be available on the Moon. Such a record would clarify not 
only the geological history of Earth, but also its chemical and 
biological history—especially since these impactors were poten- 
tially major sources of biotic precursors on early Earth (Pierazzo 
and Chyba 1999). 

Mars is presently the focus of attention with regard to the 
search for early signs of life outside Earth. Ironically, the Moon 
may be the better place to search for the remains of both early 
martian and early terran life. Most significantly, the Moon lacks 
the water capable of carrying contaminants into the interior of 
rocks through cracks. While gardening from micrometeorites is 
less severe on the surface of Mars, the deeply buried regolith on 
the Moon provides some protection. Finally, the Moon is also a 
perfect testbed for targeted sample return. 

For these reasons, we determined the likelihood that early re- 
mains of Earth, Mars, and Venus have been preserved on the 
Moon in high enough concentrations to motivate a search mis- 
sion. While others (Gladman 1997, Halliday et al. 1989) pro- 
vided estimates for the transfer efficiency and total number of 
martian meteorites impacting the Earth, there are no estimates 
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for the abundance of terran, martian, and venusian meteorites on 
the Moon. In Section 2, we consider the transfer of impact ejecta 
from Earth to the Moon’s surface immediately following an im- 
pact event. By considering both the slow moving, Earth-bound 
ejecta and the high velocity ejecta that achieve orbit around the 
Sun, we estimate the total transfer efficiency of the material. 
Through a separate numerical simulation, we estimate a rough 
transfer efficiency for venusian material. Additionally, we com- 
pute the range of expected impact velocities on the surface of 
the Moon. In Section 3 the results of Section 2 will be combined 
with mass flux estimates to calculate concentrations of terran 
and venusian meteorites on the Moon. In addition, the contri- 
bution from martian meteorites is estimated from the literature. 
Section 4 includes a discussion of the survivability of the mate- 
rial delivered to the Moon via impacts over a range of velocities. 


II. EJECTA TRANSFER FROM EARTH TO THE MOON 


Large impacts on the Earth generate ejecta with launch ve- 
locities near to or exceeding the escape velocity. This mate- 
rial interacts with the Moon in a number of ways. Ejecta with 
speeds near the escape velocity eventually fall back to the Earth, 
get swept up by the Moon, or achieve a stable orbit around the 
Earth. Ejecta with velocities exceeding the escape velocity either 
hit the Moon while escaping the Earth-Moon system or achieve 
a stable solar orbit. We explore the velocity distribution of the 
ejecta from low speeds that just allow particles to reach the lunar 
orbit up to speeds exceeding the escape velocity. We breakdown 
the velocity distribution into three regimes: 


1. Ejecta with launch speeds less than the escape velocity, but 
still large enough to reach the lunar orbit. The minimum velocity 
is determined by the Earth-Moon distance, as discussed below. 
These ejecta have a sufficiently low relative velocity with respect 
to the Moon that gravitational focusing is important. We refer 
to this as “direct transfer” and use geometrical and analytical 
methods to derive the amount of material that reaches the Moon 
via this process. 

2. Ejecta that leave the Earth-Moon system without interact- 
ing with the Moon and achieve a solar orbit, with the potential 
of interacting with the Earth and Moon at a later time. This in- 
cludes all ejecta with velocities greater than the escape velocity. 
We refer to this regime as “orbital transfer’ and use numerical 
simulations similar to those in Gladman (1997) to compute the 
transfer efficiencies. 

3. Ejecta that leave the Earth with speeds greater than the 
escape velocity and just happen to hit the Moon as they leave the 
system. These high velocity “lucky shots” are estimated along 
with the numerical calculations. 


A. Analytical Treatment of Direct Transfer 


The analysis for direct transfer treats the ejecta as a shell of 
material emanating from the surface of the Earth. This shell is 
defined by its leading and trailing edges, and the total volume 
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increases with time as the shell rises ballistically from the Earth’s 
surface. As this volume increases, the mean density of material 
within the shell is reduced. When the leading edge reaches the 
lunar orbit, the shell interacts with the Moon until the trailing 
edge also reaches the lunar orbit. Ideally, this entire process 
then repeats in reverse, ignoring the perturbing influence of the 
Moon, as the shell free falls back to the surface. 

It is important to realize the role of the lunar period in this 
transfer method. The transfer of material will be most efficient 
when the Moon is closest to the Earth, at 3.8—3.9 Gyr. The density 
of material in the shell is determined by the lunar orbit, which 
was much smaller in earlier epochs. In addition, the amount of 
mass incident on the Earth during the Late Heavy Bombardment 
(and thus the mass of ejecta thrown up by these impacts) was 
much higher. These two effects have the potential to transfer 
large amounts of terran ejecta, as described in Section 3. 

The Moon’s present mean distance from the Earth, in units of 
Earth radii, Rg, is ~60 Ro, and it is receding at 3.82 cm per year 
(Dickey et al. 1994). Paleontology gives us constraints on this 
rate as far back as 2.5 Gyr (Walker and Zahnle 1986). However, 
the rapid lunar recession in the first ~700 Myr prevents an accu- 
rate extrapolation to the Late Heavy Bombardment epoch from 
present and paleontological data. A recent study (Zharkov 2000) 
attempted to estimate an early reference point for lunar reces- 
sion. By assuming that the Moon’s present shape was “frozen- 
in” during the Late Heavy Bombardment epoch, Zharkov used 
the current data from the lunar gravity field and a model of its 
tidal interaction with the Earth to estimate this reference point. 
Zharkov estimated that the Moon’s mean distance was 21.6Rg 
about 3.9 Gyr. We will adopt this value in our calculations, and 
as the starting point for discussions of the time evolution of the 
Moon’s orbit subsequent to 3.9 Gyr. 

The leading edge of the shell is fixed by the escape velocity 
of the Earth (we consider material with velocities greater than 
the escape speed in Section 2.2). The trailing edge, however, is 
determined by the Earth—Moon distance. Thus, the maximum 
velocity is given by Vmax = 11.2 km s7! and the minimum ve- 
locity, Vmin, is determined by considering conservation of energy 
in the Earth’s potential 


1\2 
Umin = tne( 1 = z) > (1) 
n 


where vesc is the Earth’s escape velocity and n = Phon is the 
distance between the Earth and the Moon divided by the Earth’s 
radius. For the case of n = 21.6 at 3.9 Gyr, Vmin = 10.94 km s7}. 

The position of the shell is determined by solving the system 
of equations governed by the Earth’s gravitational potential, 


= ; =v, (2) 


where G is the gravitational constant, Mg is the Earth’s mass, 
and r and v are the instantaneous position and velocity of the 
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shell’s edge. This set of equations is solved separately for the 
leading and trailing edges of the shell. From this, we can deduce 
the volume of the hemispherical shell, 


E = ee (3) 


Vshe 5 
shell 3 


where Finner is the distance from the Earth to the trailing edge, and 
Touter 18 the distance to the leading edge. The density of material 
within the shell, shet, is simply the mass ejected with speeds 
between Vmin and Vmax divided by this volume. 

Chyba et al. (1994) calculated the amount of material leaving 
the surface with velocity greater than a given velocity, v, as 


M.(>v) = 0.11(p/p;)°7(v;/v)' 72m, (4) 


where p and p; are the density of the impactor and target, and v; 
is the impact velocity. In the Earth’s case the incident material is 
primarily from the asteroid belt and v ~ 14km s~! (Chyba et al. 
1994, Bottke et al. 1994). Assuming p = p; = 2860 kg m-? for 
basalt, we derive that M.(>Umin) — Me(>Umax) = 0.004m with 
n = 21.6 at 3.9 Gyr. This amount will be reduced as the Moon 
recedes from the Earth over time. 

The Moon will intersect the leading edge at a 45° angle, with 
a relative velocity, v,e, determined by the quadrature sum of 
the leading edge velocity and the lunar orbital velocity. The 
incremental path length, dP, traversed by the Moon in a time 
aT is 


dP= V1 AT, (5) 


while the incremental mass swept up by the Moon is 


2 
dM = Psheuo d P; o = An( + (==) l (6) 
v 


where A,, is the cross-sectional area of the Moon, and Vesc,m 
is the lunar escape velocity. The gravitationally enhanced cross 
section, o, is then integrated over the normalized velocity dis- 
tribution, and the incremental mass, dM, is integrated from the 
time the leading edge reaches the Moon to the time the trailing 
edge reaches it. We multiply this total mass by a factor of 2 to 
take into account the reverse process of the shell falling back 
to the Earth. Table I shows the results of the model for nine 
positions of the Earth and Moon from 3.9 Gyr to the present, 
corresponding to an Earth—Moon distance from 21.6 to 60.3 Rg. 
In particular, we report Fyirecr, the fraction of mass reaching the 
Moon, as a function of the Earth-Moon distance, in terms of the 
impactor mass. 


B. Numerical Simulations of Orbital Transfer 


A second method to transfer material from the Earth to the 
Moon is orbital transfer, which is much less sensitive to the 
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TABLE I 


The Fraction of Impactor Mass Delivered to the Moon for a Range 
of Earth-Moon Distances 


Orbital period, Interaction time, 

n Vmin» km p“ days days Fäirect 
21.6 10.94 5.9 0.30 2.0 x 1076 
26.5 10.99 8.0 0.40 1.2 x 1076 
31.3 11.02 10.2 0.50 7.9 x 1077 
36.2 11.05 12.7 0.58 5.3 x 1077 
41.0 11.06 15.3 0.66 3.9 x 1077 
45.9 11.08 18.1 0.78 2.9 x 1077 
50.7 11.09 21.1 0.84 2.2 x 1077 
55.5 11.10 24.1 0.91 1.7 x 1077 
60.3 11.11 213. 0.99 1.4x 1077 


Notes. n is the Earth—-Moon distance in Earth radii and Vmin is the velocity of 
the shell’s trailing edge. The orbital period of the Moon and its interaction time 
with the shell are given in days. Fyirecr is the fraction of the impactor mass that 
reaches the Moon. 


Earth—Moon distance. Material leaving the Earth with a veloc- 
ity greater than the escape velocity will remain in an Earth-like 
orbit as it travels around the Sun (i.e., roughly circular with 
a semi-major axis of about 1 AU). Over time, this material 
will interact with the Earth-Moon system. The goal of the nu- 
merical simulations is to determine the likelihood of impact 
with the Moon and Earth during a relatively short period of 
time. 

The dynamical studies in this paper describe the ejecta as a 
spherically symmetric distribution traveling radially away from 
the planet, similar to the assumption adopted by Gladman (1997). 
To account for a range of velocities, we simulate the particles 
with velocities, və, from 0.0 to 3.3 km s~! , which is related to 
the launch speed, v;, by 


Up = Vese + Yoo: (7) 


Since the Earth’s escape velocity is large, this corresponds to 
a relatively narrow range of velocities, from 11.2 to roughly 
11.7 km s™!. 

Each simulation consists of 225 ejecta particles, plus the nine 
planets and the Moon. The simulations were integrated using the 
dynamical integration code pkdgrav (Stadel 2001), employing a 
variable timestep to resolve close interactions with planets and 
to resolve any collisions. Since the simulation includes a large 
range of orbital timescales the maximum timestep is determined 
by the object with the shortest orbital period, which in our case 
is the Moon. To ensure detailed calculation of the particle orbits, 
the number of integrations is maximized to a time resolution of 
200 timesteps per orbit. Given that the Moon is included and 
has the shortest dynamical time (orbiting the Earth once every 
5.9 days at a distance of 21.6Rẹ), the maximum timestep is 
limited to 42 min. The minimum timestep is 1.5 min, sufficient to 
resolve any close encounters and collisions. The initial positions, 
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velocities, and masses of the nine planets were taken from the DE 
403 ephemeris, provided by the Jet Propulsion Laboratory. All of 
the simulations were run for 4775 years, each representing about 
one week of computing time on a modern desktop machine. 
The simulation results are recorded at approximately five-year 
intervals for the entire integration. 

During potential particle—planet encounters, the timestep is 
reduced to allow an accurate calculation of the interaction. The 
code uses a 16-rung multistepping ladder to reduce the timestep 
from the maximum to the minimum required to resolve the en- 
counter. Collisions are resolved at a timestep, n, by linearly ex- 
trapolating to possible encounters at the next timestep, n + 1, 
using velocities and positions of the particles at timestep n. 
Therefore, collisions with a small particle and a planet are well 
resolved, given a small enough timestep. However, errors in im- 
pact velocities can be significant. The velocity error is equal to 
ôT a, where ôT is the timestep, and a is the acceleration. At 
1 AU from the Sun, in the absence of other potentials, this ve- 
locity error is only 1.5 m s~!. However, at the Earth’s surface, it 
can be as great at 800 m s~!. To test the robustness of the code, 
we computed one simulation with half the timestep to make sure 
we recorded the same number of collisions. 

For the purposes of this study, we are interested in orbital 
transfer of material over a relatively short period of time, less 
than 5000 years. We chose this interval to bias our results toward 
material that would be rapidly buried on the lunar surface, thus 
protecting it from degradation by UV radiation, cosmic rays, 
and micrometeorite weathering. However, since the material will 
continue to be transferred even beyond 5000 years, our results 
must be considered a lower limit on the total amount of terran 
material delivered to the Moon by this process. 

Looking more closely at the reaccretion of material onto 
the Earth, Fig. 1 shows the fate of 138 impacts in the vœ = 
0.0 kms~! case, depicting the cumulative number of impacts 
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FIG. 1. The cumulative number of impacts as a function of time for the 


simulation with væ = 0.0 km s~!. Most of the impacts (123 out of 138) occur 
within the first 100 years, and the rate tapers off substantially after 1000 years. 
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TABLE II 
Impact Results for Orbital Transfer 


Scaled 
Voo km s7! Earth impacts Lunar impacts oa fo tn lunar impacts Fe Frb 
0.0 138 1 140.0 0.61 0.004 0.98 4x 1073 3x 1077 
1.0 14 0 130.4 0.06 0.0 0.11 5x 1074 5x 1077 
1.8 5 1 113.3 0.02 0.004 0.04 2x 1074 2x 1077 
2.3 1 0 101.5 0.004 0.0 0.01 4x 1075 2x 1078 
2.7 1 0 92.4 0.004 0.0 0.01 4x 1075 4x 1078 
3.3 1 0 80.0 0.004 0.0 0.01 4x 1075 5x 1077 


Notes. X® is the ratio of terran to lunar impacts, from Eq. (8) fẹ is the fraction of particles returning to the Earth, and fm is the 


N, 


m 
fraction of particles returning to the Moon. Also listed are the expected number of lunar impacts from the scaling arguments, the lunar 
transfer efficiency, F,, and the fraction of the impactor mass transferred to the Moon via this method, Fọrp. The væ = 3.3 km s7! case 


includes the ‘Lucky Shots’. 


as a function of time. We see that 89% of reaccretion collisions 
occur within the first 100 years of the simulation, and 95% of the 
collisions have occurred within 1000 years. The remaining 5% 
of the reaccretion occurs between 1000 and 5000 years. We see 
from Table II that the quantity of ejecta returning to the Earth is 
greatest for væ = 0.0 km s7! and rapidly decreases for increas- 
ing Væ. Over such a short period of time, and with simulations 
containing such small numbers of particles, lunar interactions 
with the ejecta are rare. We therefore estimate the amount of 
material interacting with the Moon from simple scaling argu- 
ments. We expect the ratio of the number of Earth impacts, Ne, 
to the number of lunar impacts, Nyoon to be 


2 2 2 
No Ro Uo + Vesc,® 
N = 2 2 2 3 (8) 
Moon Roon US0 + Vesc,m T Vpot 


where vesc, and Vesc.m are the escape velocities of the Earth and 
Moon, and Upor = Vese,g/./7 represents the contribution of the 
Moon’s proximity to the Earth’s potential well. For the distance 
to the Moon chosen for our study, 21.6Rẹ, Upo is approximately 
equal to the Moon’s escape velocity and gives a derived ratio of 
140 for væ = 0.0kms~!. This matches well with the simulation 
for væ = 0.0 kms~! and allows us to scale the other results 
accordingly. Table II indicates these scaled results, along with 
derived lunar impact fractions and transfer efficiencies, Fe. We 
use Eq. (4) to determine the amount of material launched within 
a given velocity range and use these values to derive the fraction 
of impactor mass transferred to the Moon, Fop. These will be 
used in the following section to derive the mass flux from Earth 
to the Moon. 


C. Lucky Shots 


Finally, we take into account the “lucky shots” that impact 
the Moon on their way out of the Earth-Moon system. To first 
order, for fast moving particles, the fraction hitting the Moon is 
just the cross-sectional area of the Moon divided by the surface 


area of an imaginary sphere encompassing the lunar orbit: 


(9) 


For the same conditions as the numerical simulation, namely 
Dmoon = 21.6Re, Flucky = 4.0 x 1075, the same as the transfer 
efficiency for the v = 3.3 km s7! case. Given this similarity, 
when these high velocity lucky shots are included in the mass 
flux calculation, they are treated with the high velocity particles 
from the numerical simulations. 


D. Lunar Impact Velocities 


The sections above determine the fraction of ejecta reach- 
ing the Moon. However, the velocity distribution of this mate- 
rial must be calculated to determine the likely impact velocities 
on the Moon expected during lunar encounters. Since lightly 
shocked millimeter-sized chondritic meteorites are found in the 
lunar soil (McSween 1976, Melosh 1989), there is some hope 
of intact sample survival of even fast moving impactors. 

To explore the impact velocities, we first address the veloc- 
ity distribution of the ejecta launched from the surface of the 
Earth. Equation (4) above, from Chyba et al. (1994) gives the 
integrated probability distribution in terms of the impactor mass. 
We start by taking the derivative of this equation to get the ve- 
locity distribution function for velocities within a range from v 
to v + dv, 


Py=C tals (10) 
where C is the normalization constant given by 
v= a (11) 
J omn oy aU 


We take Vmin ANd Vmax from the full range of launch velocities, 
from 10.94 to 11.7 kms~!. Thus, the probability of finding a 
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launch velocity within the range of v and v + ôv is f sa Py dv. 
Through conservation of energy in the Earth’s potential and 
knowledge of the Moon’s orbital velocity, vo-p, we can trans- 
late this into the interaction velocity at the Moon. Added to this 
is the Moon’s escape velocity, Vesc,m, to get the maximum impact 
velocity for a given launch velocity 


Nis 


Vi max = oe F vZ F Vem) , (12) 
where Vint is determined by 
1 
1 z 
Vint = (teol — ~)) A (13) 
i n 


The maximum impact velocity tells us about the total kinetic 
energy carried by the particle as it impacts the Moon’s surface. 
This corresponds to a velocity of around 5 km s™! with n = 21.6. 
However, the impactor will disperse some of this energy during 
oblique impacts, and therefore the maximum impact velocity 
does not necessarily determine the amount of shock delivered to 
the impactor. Accordingly, we compute the vertical component 
of the velocity as the upper bound on the stress induced in the 
impactor. Measuring the impact angle, 0, from the ground (i.e., 
90° is a direct impact) the vertical component is given by 
Vert = Vi max sin 0. (14) 

The probability of a given particle having an impact angle 
between 0 and 6+ 66 is given by es Po d0. Pierazzo and 
Melosh (2000a) showed that this normalized distribution is 
Po = 2sin@ cos 0 (15) 

regardless of the planet’s gravity field. Using this information, 
along with our velocity distribution, we compute the expected 


value of the vertical impact velocity, vey», over our range of 
velocities and for impact angles less than or equal to 0 


Umax fO 
Vexp(O) = T f Po Py Vvyert do dv. (16) 
v 0 


min 


Table III gives vex for a range of impact angles, and the corres- 
ponding probabilities for such angles, determined from Eq. (15). 
We can use the results from the orbital simulations to spot 
check the lunar impact calculations. We have one impact each 
for the væ = 0.0 km s7! and væ = 1.8 kms7!: one at 3.9 km s7! 
and one at 4.2 km s~! occurring 140 years and 541 min af- 
ter launch, respectively. (This last impact represents one of the 
lucky shots impacting the Moon as it leaves the system.) These 
impact velocities correspond to the respective Vimax for each 
case. Both of the impacts occurred at roughly 30° impact angles, 
corresponding to a vertical velocity of 2.0 and 2.1 km s™!, well 
within our most probable calculated vertical impact velocity. 


ARMSTRONG, WELLS, AND GONZALEZ 


TABLE III 
Lunar Impact Velocities 


Vexp km s7! 0, degrees P(<0) Impact pressure, GPa 
5.0 90 1.00 48 
2.5 77 0.95 12 
2.0 64 0.81 8 
1.0 45 0.50 2 
0.5 35 0.33 0.5 
0.1 20 0.12 0.02 


Notes. Listed are the most likely impact velocities for an angle less than 
or equal to 8, the probability that a given impact is less than or equal to 6, 
and the corresponding impact pressures experienced by the impactor on the 
Moon. 


The impact velocities also give us a rough estimate of the 
likelihood that an aggregate chunk of material will survive the 
impact in any recognizable form. Given the nature of our cal- 
culations, we want to derive the impact pressure experienced 
by the rock fragment upon impact. To get an order of magni- 
tude estimate, we use simple scaling and dimensional analysis 
to estimate this quantity. 

The total deceleration experienced by the rock fragment upon 
impact is approximately 


(17) 


where Vyen is the vertical velocity and AT = Lt is the time it 
takes to traverse the diameter, d, of the particle. We can estimate 
the mean impact pressure experienced by a particle by multiply- 
ing this deceleration by the mass of the particle and dividing by 
the cross-sectional area of the particle. Thus, we derive an order 
of magnitude estimate for the mean impact pressure which is 
independent of the size of the impactor, 


2 2 
Pave = FP Vvert> 


3 (18) 


with p = 2860 kg m~? for basalt. Table III shows the results for 
a range of likely impactor velocities. These values are in rough 
agreement with the literature-based constraints we impose for 
ejecta leaving the surface of the Earth as described in Section 4. 
It should also be noted that dissipative effects of the regolith and 
other factors may reduce the impact pressure. Thus, for the range 
of impact velocities calculated in this section, the likelihood of 
terran ejecta surviving in some large aggregate sample is quite 
high, as will be discussed in detail later. 


IHI. MASS TRANSFER RATES 


A. Mass Transfer from the Earth to the Moon 


Armed with transfer efficiencies from our analytical calcula- 
tions and numerical simulations, we can compute the mass flux 
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of terran ejecta incident on the Moon. We will proceed in four 
steps: 


1. Calculate the mass of material incident on the Earth during 
the period of interest, from 3.9 Gyr to the present, scaled from 
the lunar impact record. 

2. Determine the fraction of material that leaves the Earth 
with a given velocity during an impact event. 

3. Apply the transfer efficiencies determined in the previous 
section to estimate the mass of material that reaches the Moon. 
Ejecta with velocities between 10.9 and 11.2 kms~! correspond 
to direct transfer efficiencies. Ejecta with larger velocities cor- 
respond to either the orbital transfer or lucky shot efficiencies. 

4. Determine the fractional volume of this material in the 
lunar regolith by considering the other material accreted by the 
Moon at the same time, namely micrometeorite flux and ejecta 
from subsequent lunar impacts. 


To gain an estimate of the mass flux incident on the Earth as a 
function of time, we follow the method described in Chyba et al. 
(1994), making use of the lunar cratering record. The cumulative 
crater density for craters greater than a given size, D, due to 
impactors incident on the Moon as a function of time can be 
modeled by an equation of the form 


D —1.8 
N(t, >D) = f(t) x (=) km? Gyr™!. (19) 


4km 


We consider two possible variations in the time dependence 
of this model. First, we explore the model employed by Chyba 
etal. (1994), which describes the impact history as being roughly 
constant for the past 3.5 Gyr, and increasing exponentially for 
earlier times, 

fapt) = a(t + pe"). (20) 
We also consider the case of the lunar cataclysm, which includes 


a period of increased impact flux about 3.9 Gyr. This is modeled 
by adding a Gaussian term, 


(t-u) 


feat) = fapt) Hye , (21) 
with u = 3.9 Gyr, the peak of the cataclysm, and Te is the 
Gaussian width of the event. Each of these equations is fit to 
the lunar data provided in the BVSP (1981). Since the lunar 
data are cumulative in both crater size and time, the data must 
be fit by the cumulative number density distribution 


tage 
N.(4km) = N(t,4km) dt km~’, 


0 


(22) 


with D = 4 km as shown in Fig. 2. Similar to Chyba et al., we 
perform a x? minimization with the coefficients as free param- 
eters, leaving the decay time, t, and the width of the cataclysm, 
Te, fixed for each case (as noted in the table, t for the case of 
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FIG. 2. Analytical fits to the BVSP dataset for two models of the time 
evolution of impactor flux. The dashed line is the lunar cataclysm, and the solid 
line is the exponential model. The plot is the cumulative crater density per square 
kilometer (in both size and time) for craters larger than 4 km shown with the 
BVSP data. The effects of the lunar cataclysm can be seen around 3.9 Gyr. 


the lunar cataclysm was modified slightly from Chyba’s model 
to improve the fit). Table IV details the parameter values and 
fit results graphically depicted in Fig. 2. The values for œ and 
B given in Table IV for the exponential case differ from those 
determined by Chyba et al. (1994), but the differences are not 
significant and can be traced to differences in the interpretation 
of data from BVSP. 

Following arguments similar to Chyba et al. (1994), we now 
calculate the mass incident on the Moon as a function of time. 
The mass of an impactor can be related to the final crater diam- 
eter, 

m(D) = 0.54r v7! 8 D? p36, (23) 
where T = 1.6 x 10° (with units of kg s~!-°’) is a constant for 
lunar values of the surface gravity and surface density (Chyba 
et al. 1994), v is the impactor velocity in m s-!, De = 11,000 
is the transition from craters with simple bowl shapes to more 
complex morphologies, and D is the final crater diameter, both 
in meters. Combining Eq. (23) with Eq. (19) gives a cumulative 
crater density for objects greater than a given mass, m, 


—b 
ne >m) = fO x ( ) km? Gyr}, (24) 


m 
m(4km) 


TABLE IV 
Values of the Fit Parameters for the Lunar Cratering Record 


2 


a B y T, Gyr u, Gyr Te, Gyr x 
2.24 x 1075 9.82 x 107!! 0.144 -= — 0.50 
2.26 x 1075 3.57x107!? 1.27x 1073) 0.129 3.9 0.07 0.99 
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where b = (1.8/3.36) = 0.54. Furthermore, we can determine 
the total mass of impactors incident on the lunar surface per 
unit area as a function of time by integrating over the mass 
distribution 


Mmin on 
M(t) = i m— dm. (25) 
om 


Mmax 


Upon integration and the assumption that Mmax >> Mmin, the 
total mass of objects incident on the surface per km? per Gyr is 
given by 


b 
M(t) = fm km)? — mDr)”, (26) 


where m(D p) is the maximum mass impactor represented by the 
largest impact basin, taken by Chyba to be the Dr = 2200 km 
South Pole—Aitken crater. 

In scaling these impactor fluxes to the Earth, we take into 
account the increased gravitational focusing of the Earth through 
the ratio of impact cross sections, 


Vese,@ \2 
co _ Re TG") (27) 
Ryton 1+ (282)? 


Uj 


OMoon 


where Rg and Rmoon are the radii of the Earth and Moon, vose.g 
and Vesc,m are the escape velocities of the Earth and Moon, and v; 
is the typical velocity of the interacting particle. If the particles 
are asteroids, with typical impact speeds of 14kms~! (Chyba 
et al. 1994, Bottke et al. 1994), this ratio is ~24. Multiplying 
M(t) by the area of the Moon and taking into account the ratio 
of impact cross sections gives the total mass of material incident 
on the Earth as a function of time. 

Finally, we consider the transfer of material from the Earth to 
the Moon. In Section 2, we computed the total fraction of ma- 
terial transferred directly from the Earth to the Moon, Fäirect(t), 
as a function of time. The mass transferred via this method is 


Mair (t) = Fäirect (t)M(t). (28) 


For orbital transfer, we determined the fraction of the impactor 
transferred from the Earth to the Moon, F,,», at a distance of 
21.6Rẹ. The mass of material ejected from the Earth at speeds 
greater than the escape velocity during these impacts is pro- 
portional to the mass of the incident material, m, determined 
from Eq. (4). Using p = p; = 2860kgm~*, we derive that 
M.(> Vesc) = 0.14 m. 

Since the Earth does the bulk of the gravitational focusing 
when interacting with Sun-orbiting particles, as the Moon re- 
cedes from the Earth, the transfer efficiency will decrease in 
proportion to the ratio of gravitational potential of the Earth at 
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a distance Rọ = 21.6 Rẹ and the instantaneous value, D(t), 


Ro 


Forb (t) = For Dy’ 


(29) 


We take the instantaneous value of D(t) to have the functional 
form described by Williams (2000), 


13 (R) 


| eens 
> Ro (30) 


2/13 
D(t) = Rol ny - | , 


where (R) = —673 Rẹ Gyr! is the average lunar recession ve- 
locity, and Tef = 3.9 Gyr is the reference time for our simula- 
tions when the Moon was 21.6Rq from the Earth. (R) is deter- 
mined by taking D(0) = 60.3 Rg, the current lunar distance. 

To determine the total mass transferred to the Moon, we mul- 
tiply the ejected material by the transfer fraction 

Mor(t) = Forv(t)M (t). (31) 
Table II shows the fraction of material launched for a given væ. 
The lucky shots are included in the last integral, since the trans- 
fer efficiencies are the same as those for væ = 3.3kms~!. The 
relative contributions from orbital and direct transfer show that 
orbital transfer dominates for the past 3.9 Gyr. At 3.9 Gyr, or- 
bital transfer accounts for 58% of the total transferred mass. The 
orbital contribution almost completely dominates at the present 
epoch. 

It should be noted that these transfer methods ignore the in- 
hibiting effects produced by the Earth’s oceans and the Earth’s 
atmosphere. However, we are considering only large impactors, 
and our condition that Mmax >> Mmin Must be satisfied. To put 
this in perspective, the mass of the impactor that formed the 
South Pole—Aitken crater is roughly 1.5 x 10!° kg. If we only 
consider impactors larger than 10 km (roughly the scale height 
of the Earth’s atmosphere, and also larger than the depth of the 
Earth’s ocean), this corresponds to a mass of 1.5 x 10° kg, a 
factor of 1074 less than may. In addition, impacts of this size 
should excavate enough atmosphere for ejected material to leave 
the Earth with minimal interaction (Melosh 1988). 


B. Burial on the Lunar Surface 


As the terran material reaches the Moon, it will be buried 
by the accretion of micrometeorites and ejecta from subsequent 
lunar impacts. Love and Brownlee (1993) measured directly 
the flux of micrometeorites incident on the Earth today to be 
Micro = 4 x 10! kg Gyr7!. In the absence of other constrain- 
ing evidence, we take this to follow the time dependence of 
Chyba et al. (1994) and scale it to the Moon by considering the 
difference in surface area and gravitational focusing between the 
Earth and Moon. 

The mass of ejecta, crudely averaged over the surface of the 
Moon, can be determined from methods described in Sleep and 
Zahnle (2001). We scale the mass of ejecta produced in a given 
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impact to the volume of material removed from the crater. Cast 
in terms of the Moon’s escape velocity and the mass of the 
impactor, the relationship is 


1.3 
a") = 18.20R on )™(0, 92" ( ) m”, (32) 


esc,m 


where v; is the impact velocity on the Moon, R moon is the radius 
of the Moon in meters, p; is the density of the lunar surface, taken 
to be 2860 kg m~’, and m is the impactor mass, in kilograms. To 
determine the total amount of ejecta produced on the Moon as a 
function of time, we integrate this over the impactor distribution, 
Eq. (25), 


Mmin an 
M.j(t) = ej(t)— ; 
IE) I Mej( Ep dm (33) 


Mmax 


Again, assuming Mmax > Mmin and b = 0.54, the total mass of 
ejecta as a function of time is 


M(t) = 40.4(R moon)? #0)? 


1.3 
v (=) m(4 km)? f(t)m(D p). (34) 


i 


Thus, the total amount of ejecta available to bury terran material 
is this amount minus the amount lost to space, determined by 
Eq. (4), with appropriate lunar parameters. 

At this point, for each of the contributions of terran sources, 
micrometeorites, and subsequent lunar ejecta (collectively de- 
noted Mcont, below), we define an equivalent thickness of mate- 
rial per square kilometer per Gyr (averaged over the surface of 
the Moon) 


Meont(t) 


H(t) = — 1, (35) 
p 


with the densities, pọ, of the various materials taken to be 
2860 kg m~? for the ejecta and terran material, and 2200 kg m~’ 
for the micrometeorites. The fractional volume of terran material 
accreted subsequent to a time, tage, is determined by integrating 
each of the contributions up to the age of the surface in question 
and taking the ratio of the volume of terran material to the other 
sources 


Vier 2: o“ Herdt 
~ flage tage 3 
h i Anicro + i : Hej dt 


(36) 
Vother 


Finally, the total depth of material accreted since a time tage is 
determined by summing the contributions of all three materials, 


lage 


total = Amicro + Hej + Aer dt. (37) 


0 
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FIG. 3. Plot of the total regolith depth accreted as a function of time. The 
lines correspond to the models as stated in Fig. 2. 


Figure 3 shows Horta aS a function of time, with most of the 
material contributing to the depth of lunar regolith deposited 
during and shortly after the Late Heavy Bombardment, from 
3.9 to 3.8 Gyr. Figure 4 shows the fractional volume, or mixing 
ratio, of terran material as a function of time. Again, the highest 
fraction occurs from 3.8 to 3.9 Gyr, when the Moon’s proximity 
to the Earth facilitated transfer. 

These calculations, however, are a lower limit on the amount 
of terran material accreted by the Moon. Sleep et al. (1989), 
and subsequently Chyba et al. (1994), pointed out that, due to 
small number statistics, the Moon undersamples the largest im- 
pactors. Chyba et al. (1994) predicted more impacts on the Earth 
and argued that an impactor of mass ~1.4 x 107” kg was prob- 
able subsequent to 4.4 Gyr. Using Eq. (4) and a simple average 
of our transfer efficiencies, 1.5 x 1076 times the mass of the 
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FIG. 4. Plot of the fractional volume of terran material as a function of 
time in the lunar regolith. The lines correspond to the models as stated in Fig. 2. 
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impactor, as much as 6 x 108 kg km~? of material will reach 
the Moon, corresponding to an equivalent basalt thickness of 
0.2 km. However, it is impossible to tell when this material may 
have arrived on the Moon, depending as it does on a very few 
large impacts occurring stochastically through time. In addition, 
these impactors were likely accreted long before the end of the 
Late Heavy Bombardment, and thus before the formation of the 
lunar maria. Therefore, we ignore this contribution in the above 
calculation, but note that this very significant amount will most 
likely be buried and possibly protected by the formation of the 
maria. However, if subsequent impacts break through the maria, 
this reservoir of Earth material may be exposed to the surface. 
According to our calculations, it would be necessary to excavate 
to a depth of at least 300 m (and most likely more) to gain access 
to this layer of terran material. Assuming a crater diameter to 
excavation depth ratio of 10 to 1 (Melosh 1989), any crater 3 km 
in size or greater should provide a window to this layer. 
Finally, the effects of gardening (the breakup and mixing of 
the lunar regolith by micrometeorites and larger impactors) on 
the vertical distribution of this material will be substantial. Lunar 
cores from the Apollo missions directly indicate vertical mixing 
of up to 1 m (Mustard 1997), and samples indicating deep mixing 
to tens of kilometers have been found (Lindstrom and Lindstrom 
1986). While most of the very deep mixing occured before the 
period of interest, we can expect that the regolith will be well 
mixed to the depths indicated in this study. Therefore, any sam- 
pling of the lunar stratigraphy will have a coarse time resolution 
at best. This mixing may facilitate easier retrieval of samples 
over a range of potential ages, as subsequent small impacts will 
excavate older material directly to the surface. The amount of 
terran material on the surface of the Moon will depend largely 
on the age of the surface that is searched. Assuming the regolith 
is well mixed, we estimate the total surface abundance of terran 
to lunar material to be 7 ppm. This corresponds to ~20,000 kg of 
terran material over a centimeter deep, 10 x 10-square-km area. 


C. Mass Transfer from Venus and Mars to the Moon 


We also explored the fate of material transferred from Venus 
and Mars to the Moon. Due to Venus’ proximity to the Sun, 
the velocity required to launch a particle from the surface with 
enough energy to reach an Earth crossing orbit is Up. = 
3.3kms~—!. We performed one simulation for Venus, tracking 
225 particles for about 5000 years. After this period of time, 
only one particle hit the Earth. Thus, subject to the errors of 
small number statistics, the Venus—Earth transfer efficiency is 
~4 x 1073. Using our scaling for the gravitational cross sections 
of the Earth and Moon further reduces this transfer efficiency, 
on average, by another factor of ~10~>, giving a Venus—Moon 
transfer efficiency of 1.2 x 1077. This is 30,000 times less than 
the Earth-Moon transfer efficiency of 4 x 107° for particles 
leaving the Earth with the escape speed, indicating that Venus’ 
rocks on the Moon will be a rare find indeed. Still, an area of 
10 x 10 square km should still yield almost 1 kg of venusian 
material, if it can be identified as such. Moreover, our estimate 
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is a strict lower limit. Melosh and Tonks (1993), in a simula- 
tion of the fate of ejecta from other planets, indicated that about 
30% of the ejecta leaving Venus with a wide range of velocities 
from væ = 0.0 to 5 km s~! hits the Earth within 12 myr. This 
results in a Venus—Moon transfer efficiency of ~3 x 1076, only 
about 1000 times less than the Earth—-Moon transfer efficiency 
for particles launched from Earth with the escape speed. This 
means up to 30 kg of venusian material could exist in the same 
10 x 10-square-km area. 

The situation for Mars’ rocks is more optimistic. Estimates 
from Gladman (1997) and Halliday et al. (1989) indicated that 
15 100-g martian rocks impact the Earth each year. Scaling this 
to the Moon, and assuming the amount of martian material scales 
with the same time dependence as the rest of the incoming debris, 
we expect 6 x 1078 of the lunar material to be martian in origin, 
corresponding to about 180 kg in the same 10 x 10-square-km 
area. 


IV. SURVIVABILITY OF BIOLOGICAL 
AND GEOCHEMICAL TRACERS 


Having calculated the mass fraction of terran material on the 
Moon, we can now estimate the forms of evidence such material 
may contain. First, we identify physical constraints for specific 
evidence types. Then, we determine the planetary ejection and 
lunar impact regimes in which those constraints are satisfied. 
Finally, we convert this into a mass fraction of terran material on 
the Moon corresponding to the given evidence type for different 
impactor velocities. 


A. Constraints on Survival of Biological 
and Geochemical Evidence 


The pressures and temperatures experienced by the terran ma- 
terial during ejection from Earth and deposition on the Moon 
will determine the survival of biological and geochemical evi- 
dence. Although the impact temperatures, pressures, and fates 
of targets have been extensively considered, little information 
is available about these conditions in the projectile during and 
subsequent to impact (Pierazzo and Melosh 2000b). Given this 
scarcity, and to avoid the complications of various equations of 
state and impact parameters, we establish general and approxi- 
mate criteria to estimate the survivability of four different types 
of evidence: isotopes, significant volatile inventories, organic 
carbon, and molecular fossils (biomarkers). We emphasize that 
these calculations are intended only as an order of magnitude ap- 
proximation for the amount of material experiencing conditions 
conducive to the possible survival of these evidence types. 


1. Isotopes 


Carbon isotope systematics have been interpreted as evidence 
for life on Earth as early as 3.85 Gyr (Mojzsis et al. 1996). 
Similarly, other isotopes could also serve as biosignatures, for 
example, nitrogen or oxygen (Blake et al. 2001). Analysis of bi- 
ological isotope fractionations requires comparison of a putative 
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biological fractionation to the background signature of the host 
rock. Effectively, this imposes the condition that the terran ma- 
terial avoid significant melting. Following Melosh (1988) and 
Stoffler et al. (1991), we set the onset of melting at a pressure 
of 70 GPa, recognizing that this value actually varies according 
to specific mineral assemblages and porosities. We assume that 
most material in a 5 km s7! impact on the Moon will not attain 
peak pressures greater than this and thus should avoid substan- 
tial melting, as indicated by Ahrens and O’ Keefe (1971), Kipp 
and Grady (1996), and our own calculations (Table II). Since 
all lunar impacts due to direct and orbital transfer occur at ve- 
locities <5 km s~!, the fraction of material remaining solid is 
determined by the material ejected from Earth at pressures less 
than 70 GPa. 


2. Significant Volatile Inventories 


Tyburcezy et al. (1986) empirically determined the percent 
volatile loss of a carbonaceous chondrite as a function of the 
projectile velocity. At velocities slightly greater than 2.0 km s7! 
and peak shock pressures <30 GPa, 50% of projectile volatiles 
were lost. We therefore limit ejection to 30 GPa and lunar impact 
velocities <2.0 km s7!. 


3. Organic Carbon 


Empirical (Tingle et al. 1992, 1998, Tyburczy et al. 1986) and 
theoretical (Pierazzo and Chyba 1999) work suggests significant 
organic carbon fractions can survive impacts. In the experiments 
of Tingle et al. (1992), only small amounts of organic carbon 
were lost for impacts up to 2 kms! and pressures up to 10 GPa. 
Although other work has suggested more significant losses (e.g., 
Peterson et al. 1997), we adopt as our constraints a maximum 
pressure during ejection of 10 GPa and lunar impact velocities 
<2.0kms~!. We also note that the very fact that organic carbon 
can survive some impact conditions implies the need for rigor- 
ous criteria to distinguish extraterrestrial and terrestrial organic 
carbon. 


4. Molecular Fossils 


To our knowledge, no information is available concerning 
survival of molecular fossils under impact conditions. In fact, 
impact conditions are likely to far exceed the nonetheless se- 
vere exigencies (temperatures ~300 C, pressures ~1 GPa) of 
prehnite-pumpellyite facies metamorphism, from which 2.7 Gyr 
biological lipids have been successfully extracted (Brocks 
et al. 1999). Given the enormous informative value of biomark- 
ers and the long-term stability of hydrocarbons (Mango 1991, 
Dutkiewicz et al. 1998), we believe it is not unreasonable to 
expect that new biomarkers, corresponding to the harsh condi- 
tions of lunar impact, could usefully be defined. Toward that 
end, we calculate the fraction of terran material ejected from 
Earth at pressures <1 GPa and impacting the Moon at velocities 
<0.5 km s7!. 


193 


B. Calculation of Survival Fractions 


Until recently, the prospect that material could escape a planet 
via a natural process was considered extremely unlikely, much 
less that the material could do so without being heavily shocked 
(Melosh 1993). Experimental (Gratz et al. 1993) and observa- 
tional evidence have forced a revision of this opinion. Most 
significantly, the SNC meteorites experienced shocks no greater 
than 45 GPa, with corresponding temperatures <600°C (Stoffler 
2000, as cited by Mastrapa et al. 2001). In fact, ALH84001 
apparently traveled from the surface of Mars to Earth without 
ever exceeding 40°C (Kirschvink et al. 1997, Valley et al. 1997, 
Weiss et al. 2000). Two mechanisms might account for these 
observations. During an oblique impact, a vapor plume jet could 
entrain and accelerate surface material to the escape velocity 
with little shocking (O’ Keefe and Ahrens 1986). However, at 
low impact angle, most or all of this material will be derived 
from the impactor, while the viability of this method at higher 
impact angles (<25°) is uncertain (Pierazzo and Melosh 2000b). 
A second and more probable method, spallation, results in lightly 
shocked material due to the interference of the stress and rar- 
efaction waves near a free surface. Our calculations are limited 
to this case. 

The mass of spalled ejecta can be estimated using the method 
described in Melosh (1985), 


atl ( vi ) l 
Me(>ve) = 1|m, (38) 


PrCL Vj 2ve 


where Me is the mass of the ejecta, ve is the ejecta velocity, Pmax 
is the maximum pressure experienced by this fraction of the 
ejecta, v; is the velocity of the impactor, p; is the target density 
(2860 kg m~? for basalt), cz is the longitudinal speed of sound 
(6000 m s~!), and m is the mass of the impactor. This equation 
differs from that of Melosh (1985) only due to a mathematical 
error in the latter (for discussion and derivation of Eq. (38), see 
the Appendix). 

Implicit to Eq. (38) is the assumption that the maximum ejecta 
velocity is one-half the impactor’s. By making this conserva- 
tive assumption, we preclude the escape of spalled ejecta for a 
most probable impactor velocity of 14 km s7!; indeed, we re- 
quire more rare impacts at velocities exceeding 22.4 km s~! for 
any spall escape. This suggests that most (but not all) terran 
material on the Moon will be heavily shocked. Such material 
would nonetheless yield significant scientific returns, for exam- 
ple, by dating and characterizing the Late Heavy Bombardment 
on early Earth and in the young Solar System. Moreover, un- 
certainties pertaining to the velocity distribution of early Earth 
impactors and to maximal ejecta velocities have been treated 
conservatively in this paper. Faster impactors may have been 
more abundant 4 Gyr than now (Sleep ef al. 1989). Similarly, 
some laboratory experiments suggest ejecta velocities can actu- 
ally approach 85% of the impactor velocity (Curran et al. 1977), 
which would enable escape for a 14 km s™! impactor. 
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TABLE V 
Ejected Material as Fractions of the Impactor Mass Calculated 
for the Different Velocity and Ejecta Pressure Regimes 


vi Me/Mi M./M; M./M; Me/Mi 

(km s~!) 1 GPa 10 GPa 30 GPa 70 GPa 
14 No escape No escape No escape No escape 
22.5 1.45 x 1075 1.45 x 1074 4.34 x 1074 1.01 x 1073 
25 3.51 x 1074 3.51 x 1073 1.05 x 107? 2.46 x 107? 
30 9.14 x 1074 9.14 x 1073 2.74 x 107? 6.40 x 107? 
40 1.78 x 1073 1.78 x 107? 5.34 x 107? 1.25 x 107! 
50 2.46 x 1073 2.46 x 107? 7.38 x 107? 1.72 x 107! 
65 3.30 x 1073 3.30 x 107? 9.89 x 107? 2.31 x 107! 


Notes. vi is the impactor velocity, Me is the ejecta mass, and M; is the impactor 
mass. 


Making these conservative assumptions we determine the per- 
cent of the total terran material on the Moon that corresponds 
to the requirements set for each of the four classes of biologi- 
cal and geological evidence. Using Eq. (38), an ejecta velocity 
of 11.2 km s7}, and the pressure and lunar impact velocity con- 
straints imposed for the four material classes defined above, the 
mass of ejecta (in units of impactor mass) is derived over a range 
of impactor velocities (Table V). We then divide this number by 
the fraction of ejected material escaping Earth determined from 
Eq. (4). Once deposited on the Moon, the material delivered at 
the end of the Heavy Bombardment is expected to be buried to 
a depth of roughly 6 m within 5 myr, making protection from 
space radiation a second-order effect. Our results, in terms of 
the percentage of terran material on the Moon, are recorded in 
Table VI. For example, if all the material on the Moon were de- 
livered by an impactor with velocity of v; = 22.5 km s7!, over a 
10 x 10-square-km search area, the amounts of material having 
endured conditions permissive to the survival of molecular fos- 
sils, organics, volatiles, and isotopic signatures from early Earth 
are 0.4, 9.2, 28, and nearly 80 kg, respectively. 


TABLE VI 
Of the Terran Material on the Moon, the Mass Fraction from 
a Given Impact That Arrived There under Conditions Conducive 
to the Survival of Molecular Fossils, Organics, Significant Volatile 
Inventories, or Isotopes 


vi Biomarkers Organics Volatiles Isotopes 
km s7! % % % % 

14 No escape No escape No escape No escape 
22.5 0.002 0.05 0.14 0.40 
25 0.04 0.99 2.96 8.53 
30 0.09 2.06 6.19 17.8 
40 0.12 2.84 8.53 24.6 
50 0.12 3.01 9.02 26.0 
65 0.12 2.94 8.83 25.4 


Note. The fractions of material with lunar impact velocities of <5 km s71, 
<2 km s7!, and <0.5 km s7! are 1.0, 0.81, and 0.33, respectively (see Table MI). 
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V. CONCLUSIONS 


In this paper, we have explored the Moon as an ideal loca- 
tion to search for remnants of the early Solar System, partic- 
ularly samples of Earth not currently available to researchers. 
The Moon’s proximity and relatively unaltered surface makes 
research missions viable, and the terran samples on the surface 
are unavailable anywhere else in the Solar System. We have 
argued that terran materials are abundant and near the surface, 
with a significant fraction retaining their geochemical and bio- 
logical signatures for detailed analysis. In addition, since the 
majority of terran samples date from the end of the Late Heavy 
Bombardment, the samples in the lunar “attic” are a unique probe 
of the early conditions on Earth, and potentially contain clues to 
the earliest forms of life. 

The terran material is delivered in one of two ways, through 
direct transfer or orbital transfer. We have found that orbital 
transfer of material is more efficient for the time since the end of 
the Late Heavy Bombardment, with direct transfer of material 
only comparable in early epochs. The amount of terran material, 
7 ppm, is sufficiently large to consider a search mission. 

Before any such mission is attempted, the current stock of 
lunar material (approximately 400 kg worth) should be searched 
for terran material. Given a concentration of 7 ppm, there should 
be roughly 3 g of Earth material in the current lunar samples. 
Since lunar fines contain more than 10 million particles per gram, 
a technique of infrared spectroscopy coupled with microscopic 
imagery could distinguish hydrated silicates (common on Earth) 
from the dry, unhydrated lunar fines. Such material could then 
be isotopically analyzed to confirm its terrestrial origin. While 
this is not likely to yield much in the way of information about 
the early Earth, it would act as a proof of concept and a baseline 
for future missions. 

Recovery and identification of terran samples from the Moon 
represent a significant challenge. Still, returning a sample from 
the Moon bearing the remains of Earth life may be orders of mag- 
nitude easier than returning, say, martian samples from Mars 
bearing the remains of martian life. In this sense, lunar mis- 
sions represent an interesting proving ground for these types of 
endeavors. The risk of contamination and relative scarcity of 
terran material makes sample return missions difficult. There- 
fore, robotic missions need to be developed capable of finding 
terran material using advanced in situ measuring devices to help 
identify samples and largely driving the analysis on the Moon. 
The existence of a facility on the Moon to recover and analyze 
samples would guard against any contamination from Earth. 
However, for a complete analysis of the material, we suggest 
the best way to conduct these studies is on-site measurements 
by human observers—in essence, a return to the Moon. 


APPENDIX: JUSTIFICATION OF EQ. (38) 


The expression for the mass of spalled ejecta as a function of ejection velocity 
is based on that described and physically justified by Melosh (1985). Briefly, 
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Melosh derived an equation for the spall thickness, zs, 


_ TD 
. PrCL Ve , 


Zs (39) 


where T is the dynamic tensile strength and D is the impactor diameter. Inte- 
gration of this equation gives the volume of spalled ejecta. 

An expression for ve can be deduced from the supplementary material to 
Melosh (1985), still available from GSA, 


3 
D 
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where R is the radial distance from the impact epicenter to the point of ejection. 
Thus, zs is inversely proportional to ve and directly proportional to the cube 
of R. If a sign error is made, such that z is mistakenly assumed proportional to 
R™?, one derives the equation given in Melosh (1985). 
Avoiding that error, the volume of spalled ejecta can be determined by eval- 
uating the integral 


(40) 


Rmax 
f 2x Rzs (R)dR, (41) 
R 


min 


where the limits of integration are defined by the maximum and minimum 
ejecta velocities being considered, corresponding to Rmin = E (ays and Rmax = 
(4 2 y3 . The maximum ejecta velocity, corresponding to Rmin, has been as- 
sumed equal to half the impactor’s velocity, as in Melosh (1985). 

Evaluating this integral and substituting T = Figg (see footnote 4 in Melosh 
1985), with £ ~ 4, give the equation used in this paper. 
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